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Abstract. In this contribution wc briefly review the reverberation mapping 
technique and its results for low and intermediate luminosity AGNs. Then we 
present a monitoring campaign of high-luminosity high-redshift quasars which 
will extend these results by two orders of magnitude, probing the broad-line 
region size and black hole (BH) mass of luminous AGN at redshift ^2-3. 

1. Introduction 

The Broad-Line Region (BLR) is one of the main components of an unobscured 
AGN (ie, those where we have an unimpaired view torwards the nucleus), and 
is charactized by very broad emission lines in the UV/optical spectra (Figure 1) 
with typical widths of several thounsands km/sec due to the proximity to the 
central BH. One of the features of this region is that the observable UV/optical 
continuum flux variations are closely related to variations of the emission-lines 
fluxes. Combining the correlation between continuum and line fluxes variations 
with the widths of the lines, the mass of the central black hole can be estimated 
through the reverberation mapping technique (sec §2). 

In these days, when it is widely accepted that all massive galaxies harbor a black 
hole in their centers, reverberation has allowed real physical comparison between 
active and dormants BHs in the local universe, but so far has not probed the 
high end of the AGN luminosity range. 

In the next section we describe briefly the basic theory and results from previous 
reverberation mapping studies leading to the motivation for a high-luminosity 
high-redshift quasars monitoring campaign that we are performing since the beg- 
gining of 2005. Our main aim is to extend these results by two orders of magni- 
tude, giving the firsts measurement of the largest BHs, and thus extending our 
knowledge of the physics of AGN and their hosts into an epoch crucial for the 
understanding of galaxy evolution. 



2. Reverberation Mapping Basics 

Reverberation mapping is a technique based on the response of the BLR gas to 
changes in the central continuum source. It is based on three basic assumptions: 
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1. The ionizing continuum originates in a single central source. 

This holds because the BLR size is ~2 orders of magnitude larger than the 
continuum source. 

2. The light travel time across the BLR is the most important time scale. 

Which means that the clouds response is practically "instantaneous". In- 
deed, the recombination time is related with the electron density' hy the 
relation Tree ~ 0.1(10^° cm~^/ne) hr « 0.1 hr, which is much less than the 
typical light travel time across the BLR. 

3. There is a simple but not necessary linear relationship between the observed 
continuum and the ionizing continuum (which gives rise to the emission 
lines). 

From these points one can construct a simple ID linear model given by 

L{t) = J ^{T)C{t-T)dT (1) 

where L{t) and C{t) are the line and continuum light curves respectively. ^^'(T) 
is the transfer function (TF) which holds the information about the geometry 
and kinematics of the BLR gas. Due to technical limitations, the TF is collapsed 
to one single parameter: the time lag r between both continuum and line light 
curve variations which is obtained through a cross-correlation analysis between 
them (Peterson 2001). This lag r is used to directly measure the BLR size from 

Rblr = CT (2) 

(Peterson et al., 2004, Kaspi et al., 2005). Assuming a gravitationally bound 
system (proven to be correct for those objects with measurements from several 
lines, where the observed anticorrelation between the distance r to the BH and 
the line Doppler widths is consistent with virialized motions of the BLR gas) and 
measuring the BLR line widths a, it is possible to infer the mass of the central 
BH through 

MsH = l^^^ (3) 
where / depends on the geometry of the BLR and it is expected to be ^^1. 



3. Reverberation Results from low and intermidiate luminosity AGN 



Many important results have been found from reverberation mapping studies, 
such as the stratification of the BLR and black hole masses. Also, an empirical 
relation between the BLR size and the AGN luminosity has been determined, 
which shows that the size of the broad-line region scales with luminosity. This 
is not far from the Rblr oc L°-^ relation expected if the BLR in all AGN have 
similar rig value and ionization conditions. Kaspi et al. (2005) found that the 
scaling relation has the form 



Rblr ^0.4 



AL;^(5100 A) 
IQm erg/sec 



0.6±0.1 



pc 



(4) 
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Figure 1. Left: Optical spectrum of quasar CT260. Several UV broad- 
lines (Lya, Siiv, Civ, Cm) can be seen in the optical window indicating the 
high redshift of the object {z — 3). Right: Size-Luminosity relation for AGN 
with reverberation mapping from Kaspi et al. (2005). All our sources fall 
above the most luminous object in this diagram. 



(as it is shown in Figure lb). Preliminary results by Bentz et al (2006) suggest a 
smaller slope of about 0.52. Using -Rblr — L relation is thus possible to estimate 
the BLR radius using L alone, hence making BH mass estimation accessible 
for a large number of AGN. Such results cannot be directly applied to high-z, 
high-L sources that contain the most massive BHs, since measuring their BLR 
size requires an extrapolation by up to two orders of magnitude in luminosity in 
eq. (Q, which would be subject of a large uncertainty in the measured Rblr 
(and then in the BH mass). This is the starting point of our quasars monitoring 
program. 



4. Expanding the Luminosity Range: Quasar Monitoring Campaing 

The few attempts to extend the luminosity range of AGN with reverberation 
mapping results to high-luminosity AGNs have failed for two reasons: 1. Most 
high-L sources are very slow variables with very low amplitude (< 20%). 2. 
The line luminosity amplitude is even smaller than the continuum amplitude 
because the emission line response is averaged over the very large (~ 1 pc) BLR 
geometry. Since 2005 we are undertaking a novel strategy that uses optical 
broad band imaging to trigger spectroscopic follow up on a time scale of about 
one year. The idea is to identify high-L sources with large enough continuum 
variations that will enable accurate and meaningful measurement of emission line 
variations. Only a small fraction of high-L sources show such variations at any 
given time, hence, a large number of sources must be constantly monitored to 
identify the onset of such events. Given large enough continuum variations of a 
certain AGN, the event must be followed spectroscopically. The expected BLR 
size suggests that the spectroscopy phase should continue for at least four years. 
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Our sample consists of 56 quasars which are the brighest from the Calan-Tololo, 
2dF and SDSS catalogues spread over all RAs, with ALa(5100) > 10^^-^ erg/sec, 
redshifts that range from 2.3 to 3.3 and mean R magnitude about ~18.2. Until 
today, black hole mass estimations using reverberation mapping analysis have 
been mainly based in the variability of the H/3 (A486lA) emission line which lies 
in the optical. In our optical monitoring program we arc looking for the CiV 
(A1549A) emission line variability which, because of the the high redshifts of the 
ojects, falls into the optical window. Reverberation mapping based in the H/3 
line can be made for these objects going to the IR. 

4.1. Imaging monitoring. 

Imaging is obtained with the 0.9m Telescope at Cerro Tololo Interamerican Ob- 
servatory (CTIO) every month for a sub-sample of quasars to monitor their 
continuum emission in the R band. Monthly datasets allow us to have from 4 to 
5 light curves points in average per object in a year, enough to detect the onset 
of a high variation in their fluxes. 

We use differential photometry respect 8 to 13 field stars which allows a photom- 
etry independent of the weather conditions, since the object and the field stars 
are observed simultaneously. These flcld stars arc selected near to the quasar in 
order to appear in all the observations and have a magnitude cut-off of order of 
the quasar magnitude. To apply differential photometry we average the field star 
magnitudes in order to minimize the random variability of each of them, then 

1 ^ 

is the average field star magnitude. The rms of each field star Hght curve is also 
checked periodically to look for intrinsically variable objects. This procedure 
gives us an average magnitude almost constant within the error bars {{S) ~ 3% 
variable with a typical value of ~ 2%). The quasar Differential Light Curve 
(DLC) is given by 

.QSO \ 

S{tj) = m«^0(t,) -{S) (tj) = -2.5 log I jj== J (6) 

where S{tj) is the difference in magnitude between the quasar and the average 
of field stars at time tj (in flux it reads as the ratio of the object flux with the 
geometrical average of the field stars). So, as (S) (t) is almost flat, the time 
series of S{tj) gives us the quasar Hght curve with respect to a constant standard 
and then variations of 5 arc exclusively atributted to the quasar. To quantify 
variability we define the paremeter 

VVyar{%) = f §^ - l) • 100% (7) 

where F^ax and i^min are the maximum and minimum flux in the DLC, then 
'P'Pvar is a function of the difference Sfnin (^rnox 

(because Am oc log(-Fi/F2)). 
In this way, the Wvar parameter is just the peak-to-peak variation in flux. 
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A minimum variability is needed in the continuum to detect a correlated vari- 
ability in the emission lines (which variation is weaker because is averaged in a 
larger geometry than the continuum source). Previous reverberation campaigns 
have concluded that this minimum is ~15%, then, objects with variability larger 
than that trigger the spectroscopic follow-up. 

4.2. Spectroscopic follow-up. 

As mentioned before, we are triggering the spectroscopy follow up whenever 
we detect 15% variability using the 2.5m DuPont Telescope at Las Campanas 
Obseratory (LCO). We note that at such quasar luminosity, the typical intrinsic 
continuum variability time scale is 5-10 months. Given the (1 + z) factor, in the 
observed frame the lines are expected to follow the continuum after 1-2 years. 
During 2005-2006 we obtained spectroscopic observations for the whole sample 
to establish a zero baseline for future spectroscopy and evaluate the spectroscopic 
accuracy of our method. In 2007 we have started the actual spectroscopic follow- 
up of the first quasars. Currently, there arc 5 targets which have been followed 
up since 2007 and there are 3 more quasars reciently (2008) selected for the 
spectroscopic monitoring. 

4.3. Preliminary results. 

After 3 3'cars of monitoring we can confirm that, in general, quasars are ver)' slow 
and low amplitude variables (<0.1 mag). However, our strategy has allowed us 
to identify a subsample of 8 quasars that present a variability greater than 15% 
(about 0.2 magnitudes, shown in Figure 2) which correspond to the 14.2%. of 
the sample. The variability in some of the objects has been seen to reach up to 
20% and spectroscopic follow up has began for them. On the other hand, the 
varibility of the field stars { S ) is lower than 3%. This range is represented in 
Figure 2 using dot-dashed lines. All points of ( 5 ) fall within this region. 
Up to date DLC for our targets can be found at the website of the program 
http : //www. cec .uchile . cl/^ibotti/reverberation/. 

5. Summary and Perspectives 

We reviewed the basis of reverberation mapping technique which is a powerful 
tool to obtain black hole masses in AGN and the Size-Luminosity [Rblr — L) 
relation from which we can estimate the size of the BLR (and then estimate the 
BH mass) from a single epoch observation. Currently this relation covers almost 
5 magnitudes in luminosity, but it is subjected to large uncertainties when is 
extrapolated to high luminosity AGN. For a better determination of the slope 
of this relation we need to span the luminosity range of AGN with reverbera- 
tion mapping results of 10^'^ — 10^^ crg/'scc. To achieve that aim we present 
a monitoring campaing of a large sample of high- 2; quasars which will extend 
these results by two orders of magnitude, measuring the broad-line region size 
and black hole mass of luminous AGN at high rcdshifts. Our strategy has been 
successful identifying a subsample of high variability (>15%) quasars which are 
being followed spectroscopically since 2007. 

We expect that in the next 3-6 years reverberation mapping analysis will finally 
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Figure 2. DLCs of some of our most variable objects. Left: CTQ953 (a), 
CTQ975 (6) and CTQ320 (c). Notice tiiat the curves not only show monotonic 
flux variations, but also peaks, crucial features to apply reverberation mapping 
analysis. Right: DLCs of J224743 (a), J214355 (6) and CTQ803 (c). 



yield the highly anticipated measurement of BH masses in some of the most lu- 
minous quasars in our Universe. 
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